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We explore polaronic quantum transport in three-dimensional models of disordered organic crys- 
tals with strong coupling between electronic and vibrational degrees of freedom. By studying the 
polaron dynamics in a static disorder environment, temperature dependent mobilities are extracted 
and found to exhibit different fingerprints depending on the strength of the disorder potential. 
At low temperatures and for strong enough disorder, coherence effects induce weak localization of 
polarons. These effects are reduced with increasing temperature (thermal disorder) resulting in 
mobility increase. However at a transition temperature, phonon-assisted contributions driven by 
polaron-phonon scattering prevail, provoking a downturn of the mobility. The results provide an 
alternative scenario to discuss controversial experimental features in molecular crystals. 
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The microscopic understanding of charge transport in 
organic materials remains a tantalizing challenge after 
more than two decades of theoretical and experimen- 
tal efforts 0, Q. An established fact is however that 
transport theories developed for conventional inorganic 
semiconductors fail to properly describe organic mat- 
ter (or narrow-band systems in general) in its full com- 
plexity. One of the main reasons lies in the strong 
electron-phonon coupling which roots high technical and 
conceptual hurdles, demanding for the advancement of 
concepts and quantum transport methods able to cope 
with entangled effects of structural/chemical and ther- 
mal disorders. For strongly disordered (and soft) or- 
ganic structures (such as discotic liquid crystals, DNA 
or microcrystalline polymer semiconductors), multiscale 
approaches combining classical molecular dynamics with 
quantum transport simulations offer an interesting per- 
spective 0, H| ■ On the other hand, since the seminal 
work of Su, Schrieffer, and HeegerQ , it is recognized that 
the electron transport in 7r-conjugatcd materials encom- 
pass various complicated transport mechanisms including 
band conduction or hopping processes which can be well 
captured by the concepts of polarons. 

After early works by Holstein Q and Conwell Q , the- 
oretical progress in recent years has been achieved in the 
exploration of polaronic phenomena, including the de- 
scription of Frohlich polarons at the dielectric/organics 
interface in organic-field effect transistors (accounting for 
Coulomb interaction and screening phenomena) |8j, as 
well as the simulation of non-adiabatic polaron motionQ, 
inelastic effects on coherent polaron- like motion [l(J, or 
the revision of dynamic localization concepts in one- 
dimensional models of organic materials [111 ]. Addi- 
tionally, to investigate polaron motion in higher dimen- 
sional (3D) systems, a generalized transport methodol- 
ogy has been proposed using a mixed Holstcin-Pcierls 



model 121 in the Kubo framework. This latter work 



has provided some qualitative understanding of high 
temperature-dependent charge mobility of ultrapure or- 
ganic crystals (including anisotropic effects), but it does 
not explain the sign reversal of the bulk mobility ver- 
sus temperature derivative (^) in the lower temperature 
regime [13( which still stands as a debated issue. Indeed, 
as observed by Karl and coworkers in a-perylene sin- 
gle crystalline thin films, charge mobilities can increase 
from 1 cm 2 (Vs) _1 (10K) to around 100 cm^Vs)" 1 at 
30 K, followed by a power-law decay up to room tem- 
perature A similar mobility behavior has been also 
more recently reported in rubrene crystals-based field ef- 
fect transistors [14| and in modified pentacene polycrys- 
tallinc films with large domain sizes of over lOOnm hj. 
Notwithstanding repeated analysis, the underlying phys- 
ical mechanism drivin g su ch temperature-dependent mo- 
bility remains elusive [161 ]. Moreover, the reported clear 
evidences for band dispersion (as large as 400 meV) in 
high purity rubrene single crystals 1171 support the sce- 
nario of polaron band transport [ill, EH Eif. although the 
question of quantum interference effects [20( remains to 
be investigated and quantified in 3D disordered organic 
crystals. 

In this paper, we use a real-space order N Kubo 
methodology to theoretically explore polaron transport 
in three-dimensional disordered organic crystals, taking 
into account elastic scattering, as well as decoherence 
effects. Decoherence is introduced at a phenomeno- 
logical level, assuming that the coherent polaron mo- 
tion suffers from additional stochastic dephasing events, 
driven by thermal vibrations which play an important 
role in the low-T regime [2l|. Additionally, incoher- 
ent contributions of polaron motion (driven by polaron- 
phonon scattering [12| ) are found to dominate the high-T 
regime. The obtained temperature-dependent mobility 
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FIG. 1: (color online) (a) Typical time dependent diffusion 
coefficient for various W and temperature T. Inset: 3D crys- 
tal structure with system parameters £ mn and R m n for near- 
est neighbors, (b) Polaron mean free path versus energy for 
various disorder strengths. Only positive energies are shown 
because of symmetry. 



fingerprints are remincscent of experimental features in 
fairly ordered organic materials |13l4l5| . 

The present theoretical approach uses a polaronic 
tight-binding (TB) Hamiltonian and orthorhombic sym- 
metry (see inset of Fig. Q]), as found in rubrene, 
which is defined by H = Hp + Hw where Hp = 
0-lnO.n + Xa f^^lfix describes polarons and 



phonons. It can be obtained from a polaronic transfor- 
mation (Lang-Firsov) [ljj of a generic model for disor- 
dered polaronic materials. In thc present study, we 
use the concept of a single effective vibrational mode 
with energy uj\~12 meV, and we use as a typical di- 
mensionless electron-phonon coupling g\ = 0.7 inspired 
by available ab initio data. The effective mode used here 
is representative of the many low-frequency intermolecu- 
lar modes which are present in organic materials. The 
temperature dependent renormalized transfer integrals 
£ mn are used following the full dressing approximation 
[ijj {N\ is the Bose- Einstein function), and are given 
by e mn = £mn exp[- £ A (1 + 2N x)9\] witn £ ™n the elec- 
tronic transfer integrals (given in Fig. [T} inset). Hp gives 
rise to a gapless cosine-type temperature dependent and 
anisotropic band structure in thc 3D Brillouin zone, re- 
garded as a typical case for organic crystals [23| ■ Finally, 
Hw = Em £ m a L a m tunes the strength of an additional 
elastic (Anderson-type) disorder potential with on- 
site energies taken at random within e m S [-W/2, W/2] 
with an average modulus of e m which is W/4. Such 
uncorrelated disorder mimics impurities and defects as 
a first approximation. [l8j A refined analysis could be 
achieved by using first principles calculations and the 
precise crystal symmetry (see e.g. [24[ for details). 

The TB model is implemented into a Kubo transport 
framework which is based on the theory developed in 



Ref. [25[. This theory is generalized in the present work 
to investigate effects of static disorder. Due to the split- 
ting of the mobility into two contributions (coherent and 
incoherent) obtained in Ref. (25[, we proceed here by 
employing two complementary coherent and incoherent 
Kubo transport methodologies to explore the whole tem- 
perature range. The first one allows a real-space descrip- 
tion of polaron propagation in the coherent regime in 
the superimposed disordered (static) potential using the 
Kubo-Greenwood approach [26|. By following the po- 
laron dynamics, the elastic polaron mean free path as 
well as the weak-localization corrections to the conduc- 
tivity in the low-temperature regime can be computed. 
Second, the incoherent (phonon-assisted) contribution of 
polaron transport is conveyed by polaron-phonon scat- 
tering. It is calculated based on the theory developed 
in Ref (25|, which is here generalized to include static 
disorder effects (details below). Both contributions (co- 
herent -(-incoherent) are computed separately and add to 
evaluate the total polaron mobility. We note that there 
is no further parameter to tune the relative amplitudes 
of both. They are fixed by material parameters and the 
impact of disorder. 

From the study of the coherent polaron wave-packet 
dynamics, the elastic mean free path £ e (E) can be first 
derived from the time dependence of the diffusion coeffi- 
cient D(E,t) which writes 



t 
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taking x as the position operator along the transport 
direction, S(E — H) the spectral measure operator the 
trace of which gives the total density of states, and 

tyj) a set of random-phase states (for details see Ref. 

26j). In this approach D(t) is derived from to the 
time evolution of ((x(t) — x(0)) 2 ) driven by the opera- 
tor U(t) = 11^ exp(iHAt/h) with At the chosen time 
step. The calculations are performed for total elapsed 
computational times of i max = 61 ps. The system size 
is at least (0.24 x 0.22 x 0.14)/im 3 and periodic bound- 
ary conditions are applied. Based on this, the coher- 
ent part of the Kubo conductivity is computed using 
a(E,t = N t At) = elp(E)D(E,t)/2, with p(E) the total 
density of states. We further derive the carrier mobility 
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with fpjj (E) being the Fermi-Dirac function, c the total 
charge density fixed to 10~ 3 , and il the volume of the 
unit cell. The chosen time t is fixed by the temperature 
dependent decoherence time t v . 

By introducing static disorder, polaron wavepackets 
are elastically scattered and suffer from multiple scatter- 
ing events and quantum interferences. Elastic scattering 
first results in a diffusive regime, characterized by the 
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FIG. 2: (color online) Coherent part of the carrier mobility for 
several disorder strengths W and using a decoherence model 
with a = 2 (see text). Inset: same information using a = 3. 



polaron elastic mean free path £ e = D m&x /2v (£) max is 
the maximum diffusion coefficient and v the polaron ve- 
locity) . Fig. [T] (b) shows £ e along the transport direction, 
which follows an approximate £ e oc W~ 2 behavior for 
small W (in agreement with the Fermi golden rule) . The 
obtained values for £ e are well below the total system 
size. When £ e becomes close to the lattice constant, it 
however ceases to be a suitable characteristic length as 
the system proceeds to the strong localization regime. In 
that limit the localization length £ might be used to fur- 
ther characterize localized states. As seen in Fig. 1(b), the 
case of W = 160 meV may be regarded (in our model) 
as a transition point. One also notes that the polaron 
bandwidth varies with temperature (renormalization of 
transfer integrals). As a result, for a given static disor- 
der potential, backscattering efficiency will increase with 
temperature. 

Beyond the diffusive regime, quantum interference ef- 
fects (QIE) produce an increase of the resistance (weak 
localization). Strong (coherent) localization effects due 
to disorder in organic crystals have been studied in one 
and two-dimensional models [3| , for which all states are 
localized irrespective to the disorder strength. In three- 
dimensional models, the validity of band type conduction 
is known to be more robust with the existence of a metal- 
insulator transition, that can be tuned by the strength 
of elastic disorder |20j ■ 

In general, quantum coherence (and localization ef- 
fects) will be reduced by decoherence mechanisms pol ]. 
Decoherence mechanisms for polaron transport could be 
driven by thermal disorders including phonon-phonon, 
electron-electron as well as electron-phonon scattering 
processes. This thermal disorder includes stochastic scat- 
tering events which relate the selected coherent polaron 



with other phonon degrees of freedom. We introduce phe- 
nomenologically a coherence time which tunes the weak 
localization correction, and decreases with temperature 
following a power law. Although there is no available mi- 
croscopic theory of decoherence in inorganic (or organic) 
materials, a power-law dependence of the decoherence 
time t v oc l/T a has been derived, with a only depend- 
ing on the transport dimensionality for electron-electron 
interaction induced dephasing [20| (for temperatures be- 
low IK), in contrast to the case of electron-phonon driven 
decoherence effects for which a ~ 2 — 4 (material depen- 
dent) [2l| . We here assume such a power law for t v to be 
generic, but we take into account some variability in the 
T-dcpendent tuning of QIE by considering different val- 
ues for a. We also assume that, at 5K, t v corresponds to 
the maximum computed time {t v = i max ) and introduce 
a lower limit for t v which ensures a seamless transition 
to the semiclassical limit at elevated T. 

Fig. [J shows // coh )(T) for various disorder strengths 
W. For values as low as 40 meV (corresponding to 
elastic mean free paths in the order of 6-f nm, Fig. [1] 
(b)), no weak localization is observed within the reach 
of our maximum computed time. This can be appreci- 
ated by the absence of a time-dependent decay of the 
diffusion coefficient (sec Fig. Q] (a)). As a result, the 
low-temperature ^( coh )(T) is seen to saturate to its semi- 
classical value on the order of 200 cm 2 (Vs) _1 for the 
chosen parameters. In contrast, for a static disorder 
strength such as V^=80 meV and W^WO meV (corre- 
sponding to £ e ~ 2nm and £ e ~ 6A, respectively) sub- 
stantial contributions of QIE are obtained, as witnessed 
by the time-dependent decay of the diffusion coefficient 
(Fig. 1(a)), and result in a marked decrease of //( coh ) at 
the lowest temperatures. By increasing W by a factor of 
2 (or 4), ^( coh )(T) is reduced by 1 (or 3) orders of mag- 
nitude. The temperature-dependent decoherence effects 
then yield a positive ^ (Fig. [2]) which is independent of 
the exact details of the decoherence strength introduced 
phenomenologically. Such a behaviour can only be ob- 
served if the temperature-induced reduction of QIE dom- 
inates the band narrowing effect which is not observed for 
the negative |y found at higher temperatures in Fig. [2] 
Importantly, in contrast to common belief that a mobil- 
ity increase with temperature implies a phonon assisted 
transport regime, we obtain a positive ^ for coherent 
motion. This behavior can be expected to be relatively 
insensitive to the strength of the electron-phonon cou- 
pling or the efficiency of the phonon dressing of carriers 
assumed in the polaron Hamiltonian. 

At elevated temperatures, the phonon-assistcd in- 
coherent contribution of polaron transport (driven by 
phonon-polaron scattering) may prevail over the coher- 
ent one. By generalizing Eq. (43) of Ref. [2f| (using the 
density of states p{E) of the disordered system) , we com- 
pute the incoherent part of the polaron mobility /i( mc ) = 

2sfCr £ m t ? » / dE 1 dE 2 p(E 1 )p(E 2 )f FD (Ei)[l - 
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FIG. 3: (color online) Temperature-dependence of the total 
mobility (main frame) and incoherent part (inset) for several 
disorder strengths W and decoherence effects with ot — 2. 



f F D(E 2 )} JZo dte^- E ^/ h {eM^x{t)gl]-l} in trans- 
port direction with $ A (t) = N x e luJxt + {\ + N x )e~ iuixt and 
display its temperature dependence in Fig. [3J [27| We 
observe an activation behaviour of the phonon-assisted 
transport which is followed by a saturation regime at 
about 200K and a weak decrease up to room temperature. 
Although reminescent of the experimental data reported 
recently we discuss below why phonon-assisted 

transport is unlikely to explain the mobility downturn 
at low T. Wc plot in Fig. [3] (main frame) the to- 



tal carrier mobility ^ coil >(T) + // mc )(T), which is the 
measurable quantity in experiments. For all temper- 
atures one observes a disorder dependence of the mo- 
bility which is weaker for the high-T regime. How- 
ever, its origin in this regime is very different from the 
strong quantum interferences described in Fig. [2] at low- 
T. Importantly, we find that even for strong disorder 
the weak localization regime of coherent motion dom- 
inates the low-temperature regime of the total carrier 
mobility, resulting in ^ > for strong enough dis- 
order (not phonon-assisted). On the other hand the 
mobility decrease < 0) at elevated T is mostly 

driven by phonon-assisted polaron motion, i.e. increasing 
contribution of carrier-phonon scattering, in contrast to 
the temperature-activation concept commonly associated 
with phonon-assisted transport. 

Conclusion.-We have studied a quantum-based po- 
laron transport scenario and found that quantum in- 
terferences can play an important role for understand- 
ing temperature-dependent transport features. Such sce- 
nario is consistent with experimental trends reported in 
clean organic crystals [13l4l5j . and the obtained results 
are not within the reach of restricted semiclassical trans- 
port approaches. Although the absolute values of charge 



mobilities are comparable to what is measured exper- 
imentally, a systematic comparison with experimental 
data, including transport anisotropy, would demand the 
use of first principles calculations for the material pa- 
rameters. This work is supported by a Marie Curie Intra 
European Fellowship within the 7th European Commu- 
nity Framework Programme. 
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